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Abstract 


Periodic  structures  were  designed,  fabricated  and  characterized  to  control  the  emission  of 
electromagnetic  radiation.  These  electromagnetic  crystals  were  fabricated  using 
interferometric  lithography,  a  technique  that  lends  itself  to  large  area  periodic  structures. 
The  characterization  was  done  using  a  Fourier  Transform  Infra-red  Spectrometer. 
Extensive  rigorous  modeling  was  developed  based  on  rigorous  coupled-wave  analysis 
and  was  shown  to  provide  a  route  to  “design-to-performance”  for  these  structures. 


Electromagnetic  crystals  are  periodic  structures  that  can  be  used  to  control  the  emission 
of  electromagnetic  radiation.  The  goal  of  the  HIDE  program  is  to  directly  control  the 
emissivity  of  an  object.  The  fabrication  of  the  crystals  is  based  on  interferometric 
lithography  (see  attachment  for  a  description  of  this  method).  This  is  a  technique 
whereby  two  laser  beams  are  interfered  to  establish  a  periodic  pattern  on  a  sample.  Using 
standard  lithography  techniques,  the  pattern  is  transferred  to  the  chosen  material  system. 
Currently,  the  material  system  of  choice  is  silicon  on  sapphire  and  bulk  silicon.  In  these 
material  systems,  we  are  able  to  etch  a  periodic  pattern  of  air  holes  in  the  silicon  sample. 
In  addition,  the  material  system  can  be  metallized  to  obtain  higher  responsivity  in  a 
smaller  area.  Keep  in  mind  that  the  system  being  tested  is  simply  a  two-dimensional 
electromagnetic  crystal  that  is  finite  in  the  third  dimension.  Remarkably,  the  need  for  a 
three-dimensional  structure  to  achieve  the  desired  result  is  not  apparent.  We  have 
independently  shown  that  omni-directional  stop  bands  can  be  achieved  using  simply  a 
one-dimensional  stack  of  alternating  dielectric  material.  Although  a  “cute”  solution,  the 
utility  of  such  a  solution  remains  to  be  seen.  We  have  focused  mainly  on  a  two- 
dimensional  system. 

One  of  the  main  goals  in  the  program  was  to  “build  to  design.”  The  design  is  done  using 
rigorous  coupled  wave  analysis  (RCWA).  The  mathematical  formulation  can  be  found  in 
the  attached  paper.  The  code  has  been  verified  based  on  various  published  results  for 
slightly  metallic  and  dielectric  structures.  The  structures  that  are  being  modeled  are 
shown  in  the  figures  below.  Extending  the  technique  to  highly  metallic  structures 
produced  a  convergence  problem.  Although  qualitative  results  can  be  achieved  with 
RCWA,  quantitatively  the  results  did  not  agree  with  the  experiment.  The  experiment  was 
performed  using  a  Fourier  Transform  Infra-red  Spectrometer  (FTIR). 

Using  RCWA  as  a  design  tool,  a  variety  of  structures  were  fabricated,  characterized  and 
modeled.  Figures  1  -3  indicate  a  small  subset  of  the  data  obtained.  The  limitations  in  the 
model  become  apparent  as  the  structure  periodicity  and  aspect  ratio  are  varied.  Figures  1, 
2  and  3  shows  the  comparison  of  the  experimental  data  with  the  results  of  the  RCWA  for 
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a  periodicity  of  1.5  micron,  1.6  micron  and  1.7  micron,  respectively.  The  fill  factor  in  all 
cases  was  40%.  There  is  a  dip  in  the  transmission  response  in  both  the  calculation  and 
the  experimentally  obtained  results.  This  dip  is  due  to  the  interference  of  the  waves  as 
they  pass  through  the  sample.  By  changing  the  periodicity,  the  location  in  frequency  of 
the  dip  can  be  controlled  as  shown  in  the  comparison  between  Figure  1  and  Figure  2. 
Keep  in  mind  that  the  computations  were  done  prior  to  the  fabrication  of  the  samples, 
thus  indicating  our  ability  to  build  to  design.  Qualitatively,  the  agreement  is  very  good 
for  Figure  1 .  As  the  periodicity  of  the  structure  is  increased,  there  is  an  issue  with  the 
convergence  of  the  code.  Quantitatively  there  is  substantial  discrepancy.  The 
discrepancy  can  be  attributed  to  two  issues.  The  first  is  that  in  the  experiment  the  back 
surface  of  the  wafer  is  unpolished.  This  leads  to  a  substantial  scattering  loss  in 
transmission  due  to  the  rough  surface  on  the  back  of  the  wafer.  The  scattering  loss 
reduces  the  amount  of  energy  that  impinges  on  the  detector,  thus  indicating  a  lower 
transmission  response.  In  addition,  the  model  assumes  an  infinite  substrate,  which  will 
further  increase  this  discrepancy.  The  second  issue  is  the  convergence  of  the  code.  The 
technique  that  is  currently  being  implemented  uses  a  Fourier  basis  to  expand  the  fields  in 
the  periodic  structure.  For  metallic  structures  that  have  square  symmetry  (square  “holes” 
in  a  square  pattern),  a  few  terms  in  the  Fourier  basis  are  sufficient  to  produce  accurate 
results.  In  our  system,  we  employ  circular  symmetry  due  to  the  ease  of  fabrication.  In 
the  case  of  circular  symmetry,  a  stair-step  approximation  is  necessary  for  the  atoms.  This 
stair-step  requires  an  unrealistic  number  of  plane  waves  to  converge.  The  circular 
symmetry  requires  an  alternative  basis  to  efficiently  expand  the  fields  that  we  are 
currently  investigating.  Thus  we  observe  a  discrepancy  between  calculations  and 
experimental  results.  Note  that  as  the  periodicity  is  increased,  the  discrepancy  between 
experiment  and  computation  increases.  This  is  again  due  to  the  lack  of  convergence  in 
our  method. 

One  of  the  main  focuses  of  our  effort  has  been  to  determine  the  angular  dependence  of 
the  sample.  Figure  4  shows  the  experimental  result  of  the  azimuthal  dependence  of  the 
sample  in  reflection.  Note  that  there  is  symmetry  in  the  response  reflects  the  triangular 
symmetry  of  the  sample.  The  symmetry  is  indicated  in  the  bright  spots  that  exist  at  a 
wavelength  of  6|im.  There  are  bright  spots,  which  correspond  to  approximately  50% 
reflectivity  at  angles  of  approximately  45  ,  105  ,  165 , 225 , 285  ,  and  345  .  The  angles 
are  determined  by  the  orientation  of  the  sample  on  the  measurement  set-up.  Note  that  the 
angles  of  high  reflectivity  are  separated  by  approximately  60°.  Consider  the  wavelengths 
between  2|im  and  3|im.  The  reflectivity  of  the  sample  changes  between  40%  and  65%  as 
a  function  of  angle.  This  is  an  extremely  desirable  property.  The  reflectivity  varies  as  a 
function  of  azimuthal  angle.  Currently  we  are  investigating  controlling  the  variation  of 
the  reflectivity  as  a  function  of  the  azimuthal  angle. 
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Azimuthal  Reflectivity,  30°  incidence 


Figure  4.  Reflectivity  of  the  sample  as  a  function  of  azimuthal  angle  and  wavelength  at  an  incident  angle 
of  30°.  The  response  reflects  the  triangular  symmetry  of  the  sample. 


We  have  investigated  methods  of  improving  the  convergence  of  our  computational 
technique,  including  paralleling  the  code  to  increase  the  number  of  plane  waves  that  the 
computation  can  handle  and  identifying  various  basis  functions  that  will  more  efficiently 
reflect  the  symmetry  of  the  structure.  Fabrication  efforts  involve  investigating  multi¬ 
layer  techniques  to  implement  in  conjunction  with  the  interferometric  lithography.  These 
efforts  remain  the  subject  of  other  ongoing  investigations. 
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Gold  dots  on  Silicon 

pitch  =  1.5  microns,  diameter  =  0.6  microns 
height  =  0.1  microns 


Figure  1.  Gold  dots  on  Silicon.  The  pitch  is  designed  to  1.5  microns  and  the  diameter  of  the  holes 
is  0.6  microns.  The  height  of  the  two-dimensional  pattern  is  0.1  microns.  The  SEM  picture  is 
shown  on  the  left.  For  the  figure  on  the  right,  the  lower  curve  is  the  experimental  result  and  the 
upper  curves  are  the  calculations  for  a  variation  in  the  number  of  plane  waves  in  the  expansion. 
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Gold  dots  on  Silicon 

pitch  =  1.6  microns,  diameter  =  0.64  microns 
height  =  0.1  microns 


Figure  2.  Gold  dots  on  Silicon.  The  pitch  is  designed  to  1.6  microns  and  the  diameter  of  the  holes 
is  0.64  microns.  The  height  of  the  two-dimensional  pattern  is  0.1  microns.  The  SEM  picture  is 
shown  on  the  left.  For  the  figure  on  the  right,  the  lower  curve  is  the  experimental  result  and  the 
upper  curves  are  the  calculations  for  a  variation  in  the  number  of  plane  waves  in  the  expansion. 
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Gold  dots  on  Silicon 

pitch  =1.7  microns,  diameter  =  0.68  microns, 
height  =  0.1  microns 


Figure  3  Gold  dots  on  Silicon.  The  pitch  is  designed  to  1.7  microns  and  the  diameter  of  the  holes 
is  0.68  microns.  The  height  of  the  two-dimensional  pattern  is  0.1  microns.  The  SEM  picture  is 
shown  on  the  left.  For  the  figure  on  the  right,  the  lower  curve  is  the  experimental  result  and  the 
upper  curves  are  the  calculations  for  a  variation  in  the  number  of  plane  waves  in  the  expansion. 
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Metallic  inductive  and  capacitive  grids:  Theory  and  Experiment 


B.  K.  Minhas,  W.  Fan,  K.  Agi,  S.  R  .J.  Brueck  and  K.  J.  Malloy 
Center  for  High  Technology  Materials,  University  of  New  Mexico 
Albuquerque,  New  Mexico  87131 

ABSTRACT 

We  present  experimental  validation  of  theoretical  modeling  of  metallic  crossed  gratings 
for  the  mid  infrared  (2-5  pm)  regime.  The  gratings  are  fabricated  using  interferometric 
lithography  and  modeled  using  a  rigorous  coupled  wave  analysis.  The  electromagnetic 
properties  of  these  gratings  are  explained  in  terms  of  surface  plasmon  excitation,  which 
couples  the  incident  light  into  the  inhomogeneous  modes  of  the  gratings. 


I.  Introduction 

Crossed  gratings  (of  which  metallic  grids  are  special  cases)  find  numerous  applications  in 
the  entire  electromagnetic  spectrum.  They  are  used  in  the  microwave  regime  as  filters  and 
for  making  efficient  antennas1.  In  the  far  infrared,  they  find  applications  as  polarizers, 
beam  splitters  and  mirrors2.  Their  use  as  filters^  and  solar  selective  surfaces^  in  the 
visible  and  near  infrared  regime  has  also  been  demonstrated.  Recently  there  has  been 
interest  in  the  behavior  of  subwavelength  metallic  crossed  gratings^’2  (inductive  grids  in 
particular)  with  potential  applications  for  novel  optoelectronic  devices.  In  that  work,  the 
properties  of  metallic  inductive  grids  were  attributed  to  surface  plasmon  excitation. 
Although  this  model  explains  the  wavelength  position  of  the  experimental  features,  a 
rigorous  diffraction  model  is  necessary  to  both  accurately  describe  all  experimental 
features  and  as  a  tool  for  crossed  gratings  design. 

In  this  paper  we  present  a  rigorous  diffraction  model  validated  by  experimental  data  for 
inductive  as  well  as  capacitive  subwavelength  metallic  crossed  gratings  under  normal 
plane  wave  illumination.  We  present  the  effects  of  variation  of  grating  parameters  on  the 
mid-IR  transmission  spectrum.  To  the  best  of  our  knowledge,  this  is  the  first  rigorous 
numerical  modeling  and  experimental  verification  of  metallic  inductive  and  capacitive 
grids. 


In  the  remaining  part  of  this  letter,  we  briefly  explain  the  fabrication  and  experimental 
characterization  of  crossed  gratings.  This  is  followed  by  the  convergence  behavior  of  the 
numerical  model;  comparison  of  the  surface  plasmons  model  with  the  rigorous  diffraction 
model  and  finally  the  experimental  verification  of  the  diffraction  model. 

II.  Fabrication  and  experimental  characterization 


The  grating  used  to  validate  the  modeling  were  fabricated  using  interferometric 
lithography  and  lift-off  pattern  transfer.  A  double  polished  silicon  sample  is  used  as  a 
substrate  material,  as  it  is  transparent  in  the  (2|im  —  5pm)  regime.  After  coating  the 

substrate  with  photoresist,  interferometric  lithography^  is  used  to  expose  the  pattern.  This 
is  followed  by  developing  the  photoresist  and  metallization  by  e-beam  evaporation  of 
gold.  The  final  step  is  an  acetone  jet  lift-off  leaving  the  desired  pattern. 

Experimental  characterization  of  these  structures  is  done  using  Fourier  transform  infrared 
spectroscopy  (FTIR),  and  all  the  results  presented  in  this  paper  are  for  normal  incidence 
transmission  in  the  mid-infrared  (2  — 5  pi)  range. 


III.  Modeling 

Numerical  modeling  of  metallic  crossed  gratings  is  a  computationally  intensive  problem. 
Previously,  crossed  gratings  were  modeled  in  the  microwave  regime  by  assuming  the 
metallic  scatterers  to  have  infinite  conductivity9" While  this  approximation  works 
well  at  microwave  frequencies  where  metals  have  very  high  conductivity,  its  use  but  for 
the  visible  and  ultraviolet  is  clearly  inappropriate.  In  the  infrared,  the  suitability  of  the 
infinite  conductivity  approximation  is  uncertain.  We  use  a  rigorous  diffraction  model  to 
avoid  these  questions  at  the  cost  of  higher  computational  complexities. 

Since  the  late  1970s,  significant  work  has  been  done  in  modeling  crossed  gratings  with 
finite  conductivity.  Derrick^, 13?  Harris14  proposed  a  model  based  on  the  co-ordinate 
transformation,  Vincent15  used  a  finite  difference  method  and  Bruno16  employed  a 
variation  of  boundaries  method.  More  recently,  Kettunen  et  al.,1^  presented  numerical 
results  for  metallic  inductive  grids  based  on  the  rigorous  coupled  wave  analysis 
(RCWA)1^.  Their  modeling  was  done  principally  for  the  case  of  rectangular  scatterers 
and  included  the  dielectric  expansion  changes  suggested  by  Li1 9  to  improve  the 
convergence  of  the  algorithm.  Here,  we  present  modeling  for  finite  conductivity,  circular, 
metallic  scatterers  arranged  in  a  square  lattice  on  a  silicon  substrate.  Results  for  both 
capacitive  and  inductive  grids  are  described. 


IILa  Convergence  of  the  numerical  model 

Our  model  is  based  on  RCWA  and  incorporates  the  modifications  suggested  by  Li19. 
While  these  modifications  work  well  for  rectangular  scatterers,  in  our  case,  the 
convergence  of  the  algorithm  slows  when  the  scatterers  are  circular  metallic  patches.  We 
therefore  restrict  the  computations  to  normal  plane  wave  illumination  and  take  advantage 
of  the  symmetry  of  the  configuration^  to  reduce  the  numerical  size  of  the  problem. 
Specifically  for  normal  incidence,  the  equivalence  of  the  positive  and  negative  Fourier 
coefficients  allows  reduction  of  the  final  eigenvalue  problem  from  (2N+1)2  to  (N+l)2, 


where  N  is  the  number  of  spatial  harmonics  for  both  the  dielectric  expansion  and  the  field 
expansion.  This  allows  us  to  study  the  behavior  of  the  solution  for  reasonably  high  values 
of  N  for  an  algorithm  implemented  on  a  personal  computer. 

Validating  the  performance  of  any  numerical  algorithm  is  an  essential  part  of  establishing 
the  utility  of  the  algorithm.  Unlike  the  situation  for  one-dimensional  gratings  where  the 
convergence  properties  of  RCWA  for  metallic  gratings  are  well-established  [REFS],  little 
is  known  of  the  behavior  of  two-dimensional  RCWA  algorithms.  Initially,  we  examined 
the  numerically  calculated  diffraction  efficiency  as  a  function  of  the  mode  number  N  up 
to  the  maximum  number  of  modes  as  allowed  by  the  computational  resources.  For  the 
subwavelength  grating  of  interest  here,  the  diffraction  efficiency  is  the  ratio  of  the  power 
in  the  0th  transmitted  order  to  the  incident  power.  Figure  1  shows  the  convergence 
behavior  of  the  algorithm  for  a  capacitive  structure  on  silicon  with  the  parameters 
detailed  in  the  caption.  From  the  figure,  it  may  be  seen  that  although  the  algorithm  is  well 
behaved  its  rate  of  convergence  is  slow.  Nonetheless  the  variation  in  diffraction 
efficiency  as  a  function  of  mode  number  N  is  at  the  most  ±5%  after  N  =15.  In  the  results 


Positive  orders  along  x  or  y  direction 

Figure  1 :  Convergence  behavior  of  the  model,  A  =  1 ,20um,  d=  0.56pm,  h  =  0.1pm,  X  =  4.0  pm, 
larger  fonts,  clearer  lines,  points,  more  accurately  aligned  drawings 

shown  in  this  paper  we  have  used  N  =  15  for  circular  scatterers  and  assume  the 
convergence  error  to  behave  similarly  for  all  structures.  However,  the  final  proof  of  the 
validity  of  the  two-dimensional  RCWA  lies  in  experimental  verification.  The  remainder 

of  our  discussion  details  the  agreement  between  the  RCWA  algorithm  and  experimental 
characterization  of  actual  metallic  crossed  gratings. 


IILb  Previous  work  and  the  surface  plasmon  model 

Previous  experimental  investigations  of  metallic  sub-wavelength  crossed  grating  in  the 
infrared  have  described  their  behavior  in  terms  of  coupling  into  surface  plasmons^l.  in 
this  model,  momentum  conservation  is  used  to  derive  a  simple  design  equation^  for  the 


condition  necessary  for  normally  incident  radiation  to  couple  into  the  inhomogeneous 
modes  of  the  gratings.  This  condition  is 

AyA={emed/[(i2+j2Xem+ed)]},/2 


here  X  is  wavelength  of  incident  light,  A  is  the  pitch  of  the  grating,  i  and  j  are  integers  and 
Em  and  £d  are  the  real  part  of  the  permittivities  of  the  metal  and  dielectric  respectively. 

The  above  equation  can  be  used  to  excite  surface  plasmons  at  a  given  wavelength.  To 
illustrate  consider  the  case  of  metal-substrate  coupling  for  incident  wavelength  of  X  = 
4°|im,  this  gives  a  grating  pitch  of  A  =  1.15°|Xm  for  coupling  light  into  the  (±04:1) 
degenerate  orders  of  the  gratings.  We  verify  this  result  using  RCWA,  as  shown  in  the  Fig. 
2,  for  the  case  of  inductive  as  well  as  capacitive  grids.  Here  in  order  to  save 
computational  efforts  we  are  using  square  scatterer  instead  of  the  circular  scatterer,  as  the 
momentum  conservation  of  the  problem  does  not  depend  on  the  shape  of  the  scatterer.  As 
may  be  seen  from  figure,  the  surface  plasmons  calculations  match  very  well  with  the 
rigorous  diffraction  model.  We  also  note  launching  of  several  other  inhomogeneous 
modes  corresponding  to  different  grating  orders  as  well  as  air-metal  surface  plasmons 
coupling. 


Wavelength  (  m) 


Figure  2:  Surface  plasmons  coupling  for  the  inductive  and  capacitive  grids,  fill  factor  =  0.5|im,  h  - 
0.1  gm,  N  =  10  (a)  inductive  grids  (b)  capacitive  grids 


III.C  Capacitive  and  inductive  grids 

From  figure  2,  we  observe  attenuation  of  incident  radiation  in  case  of  capacitive  grids  and 
enhancement  of  radiation  for  the  inductive  grids.  We  experimentally  verify  this 
complementary  behavior  by  fabricating  inductive  grids  and  capacitive  grids  of  same  pitch 
and  height.  This  is  shown  in  figures  3  and  4  along  with  the  numerical  computations. 
From  the  figures  it  may  be  noted  that  a  very  good  match  is  obtained  between  rigorous 
diffraction  model  and  experiment;  furthermore,  we  observe  strong  excitation  of  surface 


plasmons  modes  for  the  inductive  grids,  as  compared  to  the  capacitive  grids.  This  is  the 
subject  of  further  study. 


Figure  3:  Inductive  grid,  air  holes  in  gold  on  Si,  A  —  1 .20(im,  d  —  0.40  |i.m,  h  —  0.1  Jim 


Wavelength  (  m) 


Figure  4:  Capacitive  grid,  gold  dots  on  Si,  A  —  1 .20(1 rn,  d  —  0.76(J.m,  h  —  0.1  |im 


Ill.d  Fill  factor  study 

Fill  factor  is  defined  as  the  fraction  of  the  area  of  unit  cell  occupied  by  the  scatterer;  for 
circular  scatterer  this  is  simply  given  as  7t/4(d/A)2.  Surface  plasmons  coupling  is  a  strong 
function  of  the  grids  fill  factor.  To  verify  this  for  capacitive  grids  we  made  samples  of 
different  fill  factors  by  keeping  the  pitch  and  height  of  the  periodic  structures  constant 
and  varying  only  the  diameter  of  the  gold  patches. 

Figures  5  shows  the  experimental  and  numerical  results  for  different  fill  factors.  It  may 
be  noted  that  the  model  and  experiment  are  much  closer  for  higher  fill  factors,  one  reason 
for  this  may  be  that  we  are  not  taking  into  account  the  finite  thickness  of  the  substrate. 
Also  from  the  figure  we  see  that  the  surface  plasmons  coupling  get  stronger  for  higher  fill 
factors  and  reaches  a  maximum  value  after  a  certain  fill  factor.  We  expect  similar 
behavior  for  the  case  of  inductive  grids. 
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Figure  5:  Fill  factor  variations,  A  =  1 .2(im,  h  -  0.1  (im  -  0.82  |im,  X  -  4.1 5  pm 


IV.  Conclusion 

We  have  presented  experimental  and  numerical  results  for  metallic  inductive  and 
capacitive  grids.  A  physical  model  based  on  surface  plasmons  coupling  is  used  to  study 
the  behavior  of  these  grids.  This  model  is  supplemented  with  a  rigorous  diffraction  model 
to  study  effects  of  various  grating  parameters.  Experimental  and  numerical  evidence 
shows  that  the  surface  plasmons  coupling  is  a  strong  function  of  the  grids  fill  factor;  it 
increases  with  bigger  fill  factor.  We  will  present  further  details  on  the  behavior  of  these 
grids  in  a  future  paper. 
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